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We suggest a duality between the standard (dynamical) and statistical distri- 
butions of partons in the nucleons. The temperature parameter entering into the 
statistical form for the quark distributions is estimated. It is found that this effective 
temperature is practically the same for the dependence on longitudinal and trans- 
verse momenta and, in turn, it is close to the freeze-out temperature in high energy 
heavy- ion collisions. 

The description of hadron production using statistical models has been pioneered several 
decades ago by E.Fermi I.Pomeranchuk [2], L.D. Landau [S] and R.Hagedorn 4|. The 
transverse momentum spectrum of particles produced in hadron-hadron collisions can be 
been presented in the simple form ~ exp(— m^/T), where —m^t is the transverse mass 
of the hadron h and T is sometimes called the thermal freeze-out temperature. 

As is well known, the statistical (thermal) models have been applied successfully to 
describe hadronic yields produced in heavy-ion collisions (see, for example, and refer- 

ences therein). The temperature obtained in these analysis is often referred to as chemical 
freeze- out temperature and is consistently slightly larger than the thermal freeze-out tem- 
perature. At the same time, the source of very fast thermalization is currently unknown 
and alternative or complementary possibilities to explain the thermal spectra are of much 
interest. 

Situations where statistical models have been applied while the notion of statistical sys- 
tem was not obvious are not uncommon. In particular, the statistical model was successfully 
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applied to analyze deep inelastic lepton-nucleon scattering (DIS) [10|, where the statisti 



cal form 



12Hl4| for distributions of unpolarized and polarized quarks and gluons (partons) 



in a nucleon was used. 

This model may be compared with more standard parametrization of parton distributions 
based on the Regge theory at low x and quark counting rules at large x [l6]- [20 1 . It is 
especially interesting to find a counterpart and give a physical meaning of the temperature 
introduced in the statistical model for the description of parton distribution functions. This 
is the main goal of this paper. Making this comparison we suggest a concept of the duality 
between the statistical and dynamical descriptions of the DIS and strong interactions of 
particles and explore its consequences and implications. We do not consider the question 
of scaling violations for different values of Q 2 , these will come into play when considering 
gluon or quark radiation and need higher order diagrams to play a role. 

The standard general form of parton distributions at the input energy scale Q%, which is 
used in the perturbative QCD |21] as the initial condition for the Q 2 evolution, describing 
the experimental data on the DIS (see, e.g., 16] ) is of the following general form: 



xq{x) 



A s x~ aa (1 - x) bs + B v x av (l 



x 



(1) 



where the first term is a generic description of the sea quark distribution which is enhanced 
at low x, according to the DIS experimental data, whereas the second term is due to the 
valence quark distribution. Note that to fit the experimental data on the DIS a much more 
complicated form for xq(x) is used in the form of PDF's [njj-j^. 

In particular, the valence part at x < l/b Qv can, to a good approximation, be written as: 



xq v (x) ~ B qv x av exp (-xb qv ) . 
Introducing the quantity x = l/b qv Eq.(j2]) can be presented also in the following form 



(2) 



xq v (x) 



B qv x av exp 



(3) 



It would be interesting to find a link of the parameter x to the temperature T. Therefore, 
let us present a short review on the modern understanding of the statistical model of the 
hadron production and the nucleon structure functions. 

A statistical approach was applied in great detail to study the quark structure of nucleons 



in Ref. 



10], using earlier approaches as a starting point [l2j, |l3j. The nucleon is viewed as a 
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gas of massless partons (quarks, antiquarks, gluons) in equilibrium at a given temperature 
in a finite size volume. According to the statistical approach, the quark distribution in a 
proton q(x) in the infinite momentum frame (IMF) at the input energy scale Qq is fitted as 
a function of the light cone variable x in the following form: 

= B^x^ B 2 X 0q x b > 

Xq[X) exp ((x - X 0g )/x) + 1 exp(x/x) + 1 ' U 

where x plays the role of a temperature inside the proton and Xo q is the chemical potential 
of the quark inside the proton. All the parameters are found from the description of the 



DIS. The same form is suggested in 11( for antiquarks xq(x) by changing the sign for the 
chemical potential Xq q . Actually, Eq.Q is not the exact statistical form, this is similar to 
the Fermi-Dirac distribution. A meaning of the factors x bl and x b2 in Eq.Q is the same as 
the meaning of factors x av and x~ aa entered into Eq.([T]). They are dictated by the Regge 
behaviour at x — > and the DIS experimental data which show some enhancement for the 
see quark distribution xq s (x) in comparison to the valence quark distribution xq v (x) at low 
x. 

Note that the statistical weight of the quarks can be written in the form exp((E q — fi)/T), 
where E q = (P ■ p q )/m is a quark energy in the nucleon rest frame, p q and P are the four- 
momenta of quark and nucleon, respectively, m is the nucleon mass. 

Passing to the infinite momentum frame the quark distribution in a nucleon over the longi- 



tudinal momentum p z can be represented in the following statistical form 13| : 



q(p z ) ~ exp { - ( Et _S PZ ~ £ ) \ » ( 5 ) 



V vnT T , 

where E, P and m are the energy, the momentum and the mass of the nucleon moving in 
the z-direction, e and p z are the energy and momentum in the ^-direction of a parton in the 
nucleon, T is the temperature and \x is the chemical potential. For massless partons one gets 



q(x) ~ exp (^-— j = exp J = ex P(") ■ ( 6 ) 

where x = p z /P is the longitudinal fraction of the parton momentum and x = 2T/m. This 

nn 

form agrees with the result obtained in [12|, [13J. The inclusion of of the quark transverse 
momentum k qt results in Ee — Pp z ~ m 2 x(l + k qt /(x 2 m 2 )). It means that the relativistic 



invariant variable x is replaced by the variable x(l + k qt /(x 2 m 2 )) [23] 
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This form is similar to Eq.Q if x = l/b v . In this case the proton temperature is T = 
x ■ m/2 = m/(26„((5o)) for current quarks which are practically massless. Approximately 
one has for the valence w-quarks b Uv ~ 2.5 — 3 and for d-quarks bd v — 3.5 — 4 at Ql = 2 — 4 
(GeV/c) 



16 



19] ; therefore, T ~ 120 — 150 MeV for the massless quarks. 



The similarity of the thermal form for the quark distribution given by Eq.([6]) and the 
dynamic form for q v (Eq.pj) maybe interpreted as a duality of the thermal and dynamical 
descriptions of the parton distribution in proton. This property is by no means surprising 
as both parametrization are fitted to describe the same data. In Fig.l, the valence w-quark 
distribution in the proton as a function of x is presented, the solid curve corresponds to the 
CKMT parametrization [la], see the second term of Eq.ffl at Ql = 4:(GeV/c) 2 ; the dotted 
line corresponds to the statistical (thermal) BS model [ill ], see the first term of Eq.(jl]) at the 
same value of the dash-dotted curve corresponds to the statistic parametrization given 
by Eq.([2l). Figure 1 shows that all three lines are very close to each other at 0.01 < x < 0.4. 

Let us compare the parameters of the statistical parametrization for u v (x) given by Eq.(]2]) 
and the BS parametrization (Eq.(jl])). Equation (00) is not a thermal distribution. The 
parametrization given by Eq.(T5]) is different . There is thus no contradiction in getting 50 
MeV with one parametrization and 120 - 150 MeV with another one. While the temperature 
has an evident dynamical counterpart T = m/(2b v ) , a similar relation for the chemical 
potential X 0q of valence quarks is not obvious. The temperature T = m/(2b v ) entering into 
Eq.([2]) is about 120 — 150 MeV, whereas the model of Ref. [ll( suggests that T is about 50 
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13|. 



MeV and the chemical potential for w-quarks Xq u ~ 216 MeV as first found in 

The quark distribution function in a proton using the variable x and the transverse 
momentum k t can be represented in the factorized form f q (x, k t ) = q(x)g(k t ) which of 



course cannot be valid at all values of x 
lattice QCD 



24j . This is supported by simulations within the 



25] and by the observation of the so-called "seagull" effect [26[ , i.e., the weak 



rc-dependence of the mean transverse momentum < p t > of hadrons produced in hadron 
inclusive reactions at low x, e.g., x < 0.5 and the strong x-dependence of < p t > (x) at 
x — > 1, see for example [^ and references therein. 

Developing the approach of the previous section one can fit the quark distribution g(k t ) 
also in the statistical form, i.e., 



exp( 



(-kt • 



(7) 
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where €kt = ykf + 171^ is the transverse energy of quarks in proton and m q is the quark 
mass. For massless quarks (the current quarks) g{k t ) can be represented in the following 
form: 

g{k t ) ~ exp(-^) = exp(- < ^ - ). (8) 

Applying here the same effective temperature T ~ 120 — 150 MeV, as for longitudinal 
momentum distribution, we get similar results on < k t >, which were obtained for valence 
current quarks (see e.g. [23J and Ref. therein). Actually, these values for < k t > of quarks 
in proton are used to get the inclusive transverse momentum s pec tra of hadrons produced 



in p — p collisions at not large p t , which were obtained in 



27b 



model (DPM) [30] or the quark-gluon string model (QGSM) [3J| 



29( within the dual parton 



Therefore, the transverse momentum distribution of partons in proton can also be de- 
scribed in the statistical form with the same value of the temperature T as the parton 
distributions in proton over the longitudinal momentum or over its fraction x. 

Note that the longitudinal and transverse momentum dependences have been related in 
a model approach implying Lorentz (rotational) invariance 23 1 
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We analyzed the parton distribution functions in a proton using the statistical model [10 , 



131 ] . We suggest a duality principle which means a similarity of the thermal distributions 



of partons (quarks and gluons) and the dynamical description of these partons. This duality 
allowed us to find an effective temperature T ~ 120 — 150 MeV for the massless quarks. 
Concluding our study of this problem one can suggest another interpretation for 
this effective temperature T entered into the quark distributions given by Eq.([6]) 
and Eq.flT]) as the effective widths for the quark distributions over the longitudinal 
and transverse momentum, the values of which approximately are the same. 

The inclusive spectra of hadrons produced in NN collisions at high energies are not 
described by the thermal statistical distribution. However, using the thermal statistical 
form for quark distribution as a function of the longitudinal momentum fraction x and the 
internal transverse momentum k t one can get a satisfactory description of the experimental 
data on these spectra. 

It happens that the freeze-out temperature in central heavy-ion collisions at zero chem- 
ical potential or the effective width of the energy distribution of pions has a similar 



value which was estimated in this paper for the massless quarks in a free nucleon 



32j-|34]. 
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FIG. 1: The valence it-quark distribution in the proton as a function of the Bjorken variable x at 
Qq = A(GeV/c) 2 . The solid curve is the CKMT parametrization [16J, the dashed line corresponds 
to the statistical parametrization give n by Eq.@ the dotted curve corresponds to the BS statistical 



(thermal) model suggested in 



10, 



11] 



FIGURE CAPTIONS 



1. The valence w-quark distribution in the proton as a function of the Bjorken variable 
x at Qq = A(GeV/c) 2 . The solid curve is the CKMT parametrization [16J, the dashed 
line corresponds to the statistical parametrization given by Eq.(H]) the dotted curve 
corresponds to the BS statistical (thermal) model suggested in (lCJ, 



